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The electronic properties of monolayer phosphorene are exotic due to its puckered structure and 
large intrinsic direct band gap. We derive and discuss its band structure in the presence of a 
perpendicular magnetic field. Further, we evaluate the magneto-optical Hall and longitudinal optical 
conductivities, as functions of temperature, magnetic field, and Fermi energy, and show that they are 
strongly influenced by the magnetic field. The imaginary part of the former and the real part of the 
latter exhibit regular interband oscillations as functions of the frequency u> in the range hw ~ 1.5 — 2 
eV. Strong intrahand responses in the latter and week ones in the former occur at much lower 
frequencies. The magneto-optical response can be tuned in the microwave-to-terahertz and visible 
frequency ranges in contrast with a conventional two-dimensional electron gas or graphene in which 
the response is limited to the terahertz regime. This ability to isolate carriers in an anisotropic 
structure may make phosphorene a promising candidate for new optical devices. 

PACS numbers: 78.20.Ls, 42.70.-a, 73.43.-f, 81.05.Zx 


I. INTRODUCTION 

Graphene possesses extraordinary properties but its 
application in device fabrication is limited by the zero 
band gap. Graphene-based transistors suffer from a low 
on-off current ratio due to its gapless structure Hi- 
This has lead to an intensive search for materials with 
a finite band gap including silicene Q, germanene 
M 0 S 2 , and other group-VI transition-metal dichalco- 
genides Q. Despite the fact that these materials show 
a high on-off ratio, their carrier mobility is considerably 
lower than that of graphene and restricts their applica¬ 
tions in electronics and optoelectronics Moreover, 
a very large intrinsic direct band gap in these mate¬ 
rials, render them unsuitable for near-infrared optical 
telecommunication and mid-infrared applications. Thus 
the search continues for a two-dimensional (2D) semi¬ 
conducting material, with a direct band gap, high carrier 
mobility, and with the potential to form excellent con¬ 
tacts with known electrode materials. 

Recent developments in the experimental realization 
of 2D phosphorene has attracted much interest due to 
its potential applications Unlike graphene, M 0 S 2 

and other related materials, electrons in phosphorene are 
highly dispersive and delocalized along the out-of-plane 
direction [I^. Phosphorene has a honeycomb structure 
of black phosphorous with large intrinsic direct band gap 
of 1.52 eV assessed by tight-binding m and density func¬ 
tional theory [la - ll^ . Gompared to graphene, the puck¬ 
ered structure of phosphorene exhibits lower symmetry 
that gives rise to in-plane anisotropic properties in mo¬ 
mentum space. Garrier mobilities are very high at room 
temperature and exhibit a strongly anisotropic behaviour 
in phosphorene-based transistors with a high on-off ra¬ 
tio [3, 0, HI]. Due to its highly dispersive band struc¬ 
ture, phosphorene exhibits a high carrier mobility and a 
large optical conductivity [l^- A linear dichroism has 
been shown in the computed absorption spectra, in that 
the positions of the lowest energy absorption peaks differ 


strongly for the two in-plane directions [T^. In addition 
to the traditional mid- and near-infrared emission, the 
unique anisotropic electronic and photonic properties of 
phosphorene may allow for the realization of novel optical 
components such as polarization sensors and anisotropic 
plasmonic devices |15l |. Already due to its high optical 
efficiency [l6l| . phosphorene has shown a high potential 
for optical device applications [i3- 

Optical transport properties of graphene show 
good agreement between theory and experiments 1181 . 
Magneto-optical properties of topological insulators [1^ 
and other single-layer materials, such as M 0 S 2 [1^ and 
silicene [2l| , have also been studied. Landau levels (LLs) 
are formed in the presence of an external magnetic field. 
Transitions between the LLs generate absorption lines in 
the magneto-optical conductivity and were used to 
distinguish between the topological insulator and normal 
(band) insulator phases in silicene in the presence of spin- 
orbit interaction and staggered potential [ll]. From a 
fundamental point of view, many efforts have been made 
to explore different properties of phosphorene at zero 
magnetic field whereas limited work has appeared for fi¬ 
nite magnetic field . Accordingly, studies of magneto¬ 
optical properties are timely and expected to increase our 
understanding of this material. As will be shown, an im¬ 
portant difference with graphene and other 2D systems 
is that their magneto-optical response occurs in the tera¬ 
hertz (THz) regime whereas phosphorene’s can be tuned 
in the microwave-to-THz and visible frequency ranges. 

In this work we study the magneto-optical transport 
properties of monolayer phosphorene. In Sec. II we de¬ 
rive and discuss its band structure in the presence of a 
perpendicular magnetic field. Further, using Kubo-type 
formulas in Sec. Ill we evaluate the optical Hall and lon¬ 
gitudinal conductivities. We proceed with a discussion of 
the results and of the power absorption spectrum in Sec. 
IV. We then show briefly the oscillator strength of the 
optical transitions in Sec. V and conclude in Sec. VI. 
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II. BASIC EXPRESSIONS 

We start with the widely used two-band model for 2D 
phosphorene [ll|) 113 1 in which the low-energy Hamilto¬ 
nian is 


/ E<^ + {a'\ll+[inl)/2 0 \ 

H=\ , ( 1 ) 

V 0 E>^-{x'nl + in^l)i2 j 

where a' = a + /Eg, A' = A -I- /Eg, 7 = 8.5 * 10® 
m/s,a = l/niex = l/1.47me, /3 = l/mf,y = l/0.83me, 
A = 1/mhx = 1/10. 66 me, rj = l/rntx = l/1.12me, and 
Eg = E^ — E^ = 1.52 eA<^. The minimum of the conduc¬ 
tion band occurs at E^ = 0.34 eV and the maximum of 
the valence band at E^ = —1.18 eV. Further, 11 = p-l-eA 
is the 2D canonical momentum with vector potential A. 
Using the Landau gauge A = (0,i3a:, 0) and diagonaliz¬ 
ing the Hamiltonian (1), we obtain the eigenvalues 


= F;'* + s{n -H 1/2)JEj/, n = 0,l, 2, 3,... (2) 

where s = -|- 1 (—1) for electrons (holes), E^ = E^/^, 
uj/ = with uj/ = eB/= 2.696ujc, ujj/ = 
eB/^m/^mhy = 2.2076wc, = 1/a', = 1/A', 

and Wc = eB/rrif,. It is interesting to note that unlike the 
anisotropic zero magnetic field dispersion, the LL spec¬ 
trum is independent of the in-plane wave vectors. The 
corresponding normalized eigenfunctions are 


'^n,ky{r) 




(3) 


where u" = — xo)/l and I = 

yJh/eB is the magnetic length and 4>n{u) are the har¬ 
monic oscillator functions. If we use the Landau gauge 
A = (—By, 0,0), the eigenvalues are again given by Eq. 
(2) and the eigenfunctions are obtained from Eq. (3) by 
replacing x with y wherever they appear. 

The density of states (DOS) is given by 

^ E - ^n), (4) 

n,s,ky 


where 5*0 is the area of the system. The sum over ky in 
Eq. (4) is evaluated using the prescription (feg = Lxl2l‘^) 
Y.ky {Ly/2T:)gs j^^dky = (S'o/Do)^^, where Dq = 
27rP and gs = 2 is the spin degeneracy. The Fermi energy 
Ep is determined from the electron concentration ric, 

D(E)f(E)dE = gjDoY,fiE:), (5) 

n,s 



where the Fermi-Dirac distribution function is written as 
/(E®) = (1 -H exp[/3(F;® - E)])-^ with /3 = l/ksT. In 
Fig. I (top panel) the magenta solid curve shows Ep 
obtained numerically from Eq. (5) as a function of B for 


a realistic value of the electron density [ 23 , ric = 1 x 10^® 
m“^, together with the LLs obtained from Eq. (2). 

Assuming a Gaussian broadening of the LLs, Eq. 
(4) is rewritten as D(E) = (gs/Dc) J2n,s [ “ (^ “ 

E®)^/2F^], where Dc = DoFy^ and F oc -x/H is the 
width of the Gaussian distribution [^. The dimension¬ 
less DOS at E = Ep, D(Ep) = D(B), is shown in Fig. 
1 (bottom panel) as a function of the magnetic field, for 
F = 0.2^/B meV (black solid curve) and F = 0.1^/B meV 
(red dotted curve). For weak fields B the level broaden¬ 
ing is important whereas for large fields it may become 
smaller due to the ^/~B dependence since the distance be¬ 
tween the LLs increases linearly with B. 



FIG. 1. Top panel: Fermi energy (solid magenta) as a function 
of the magnetic field for T = 1 K and density ric = 1 x 10^® 
m“^. The dashed curves are the LLs evaluated from Eq. (2). 
Bottom panel: Dimensionless density of states as a function of 
the magnetic field B for two different values of the LL width. 


III. LINEAR RESPONSE CONDUCTIVITIES 

We consider a many-body system described by the 
Hamiltonian H = Hq -|- Hi — R • F(t), where Hq is 
the unperturbed part, Hi is a binary-type interaction 
(e.g., between electrons and impurities or phonons), and 
—R • F(t) is the interaction of the system with the ex¬ 
ternal field F{t) [ 13 . For conductivity problems we have 
F(t) = eE(t), where E(t) is the electric field, e the elec¬ 
tron charge, R = is the position operator 

of electron i. In the representation in which Hq is diag¬ 
onal the many-body density operator p = has a 

diagonal part and a nondiagonal part For weak 
electric fields and weak scattering potentials, for which 
















3 


the first Born approximation applies, the conductivity 
tensor has a diagonal part and a nondiagonal part 
part, n,iy = x, y. 

In general we have two kinds of currents, diffusive and 
hopping, but usually only one of them is present. In 
the present problem, due to the magnetic field we have 
only the hopping current since the diffusive one vanishes 
due to = 0, see Eq. (6). For elastic scattering the 
diffusive part of is given by 


ydi/ 


(w) = 


/3e^ 




( 6 ) 


nd 
^ fiiy 


(w) = 


ihe'^ 




- Eqi){E^ - Ec_i + hio + lE) ’ 


(7) 


y,n ,n 


— Uy ^■\/n+ lSn\n+l “t“ -x/n^riyn—l) ^ky,k' 


(9) 


X. and - 


nj-i- After making the changes n—1—^-nin the 


where AJ’® = = v^. — v 

E; 

first sum, we combine the two sums and obtain 

[fu+i- 




s,s' ,n=0 


(n + I) 


L/: 


n+l,n' 


(7^+i,n + Exi + iT) 


[fn- 


where Sq = L^Ly, is the relaxation time, ui the fre¬ 
quency, and the diagonal matrix elements of the ve¬ 
locity operator. Further, = [1 -|-exp/3(£'^ — Ep)]~^ is 
the Fermi-Dirac distribution function with (3 = l/fc^T, 
T the temperature and Ep the Fermi level. 

The collisional part cr'^Jiuj) of cr^j^ is much smaller 
than the nondiagonal part and we neglect it. As 

for given in Ref. [^, one can use the identity 

/c(l - /c')[l - exp/3(£’^ - Et:;,)] = fc - fc and cast [H] 
its original form in the more familiar one 


7::;+i(C„+i+?*a; + zr)J 

In the limit F —>■ 0, w —>■ 0 and s = s' Eq. (II) 
yields zero. The matrix elements of the velocity op¬ 
erators are nonzero only for n' = nil. Regarding 
the sums over s, s' for convenience we write = 

!]+,+ + S-,- + S-r.- + S-.-f subscript +/- 

denotes the conduction/valence band. After performing 
the summation over s and s' we obtain the absorption 
spectrum of the real part of in the form 




= ttcto y^(n i 1) 


n—0 


where the sum runs over all quantum numbers IC) = 
|n, s,fcy) and |C') = |n',s',fcy) provided C ^ C■ The 
infinitesimal quantity e in the original form has been re¬ 
placed by F to account for the broadening of the energy 
levels. Equation (7) has been recently applied to phos- 
phorene in the absence of a magnetic field in Ref. Halbl . 
The evaluation of and is outlined in the Ap¬ 
pendix. Using Eqs. (1) and (3) we obtain 

'^x,n,n' ~ ^^x l,n' ~ Vn -|- 16n+l,n'') ^ky,k'y-: (8) 






[/n + l-/n]'S(j/ 


l/(A 


-)2 


huj^ / (Aj. 


[/n~-//+l]'S(Jn;n+l-rfi^) 




l/(A 




.( 12 ) 


Here and 7r(5(a;) = V/^x^+V'^). Using 

Eqs. (8) and (9) into Eq. (7), and carrying out the sum 
over n' similar to Eq. (10) and then making the changes 
n — lA-mA-nm the first sum, we obtain the optical 
Hall conductivity as 

[fu+l - 




E 

s,s' ,n—0 


(n-f 1) 


1S,S 

^xy 


tS,S' 


(7n+l,n + Ex! + ir) 


[f^ - /4i]^: 


s,s 

xy 


^n^n+li^n^n+l + ^ + ^F)-! 


(13) 


where A*’® = and A*’® = {v^-v^)*{/3m'^y^v^- 

Vrn'hxV^)- In the w = F = 0 and s = s' = e limit Eg (13) 
yields the quantized Hall conductivity of a 2DEG 


where and = (m',,w//^e)(/3/^)- The 

results for s ^ s' are given in the appendix. Since |(/) = 
|n, s,fcy), there will be one summation over ky which, 
with periodic boundary conditions for ky , gives the factor 
So/2pP. Substituting Eq. (8) into Eq. (7), summing 
over n', and setting cto = —gshe^/2pl^, we obtain the 
longitudinal nondiagonal conductivity as 


<x{^) = * 0-0 


E 




„=i +huj + iT) 


+ icTo 




s,s' ,n—0 


^n,n+l 


(7^’,n-|-l + fiw -I- ^F) 


,( 10 ) 


2 

(. 0 ) = ~ ■ 

n—0 


(14) 


Also, apart from the absence of valley degeneracy and the 
appearance of (n-|-1) instead of (n-|- 1/2), the dc limit of 
Eq. (13) gives results similar to those of graphene [^ . 
Following the procedure adopted for the component 
and using Eq. (10) we obtain 

lma2y{uj) = TTCTo '^{n + 1) 


[fi-fi+l] S{hul-hu,) [/--/-+i] S{huj^-huj) 


hto^jA] 


TTT 


huj^ jAj. 


, [/n“-//■+!] 




/A. 




y\! 


X 


.(15) 
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IV. DISCUSSION OF RESULTS 

The energies of the positive branch levels in Eq. (2) are 
different from those of the negative branch due to the dif¬ 
ference in the E® values and cyclotron frequency in each 
band (electron/hole). Due to Huj^ << E®, the intraband 
and interhand optical transitions belong to two widely 
separated regimes: the former is in the microwave-to- 
THz range and the latter in the visible frequency rairge. 
We will first consider the latter that involve transitions 
between neighbouring LLs (n' = n± 1) and s ^ s'. In all 
results shown below the common parameter are temper¬ 
ature T = 10 K and level broadening T = 0.2'/B meV. 

0.360 
0.355 
> 0.350 

lU 

uq 0.345 
0.340 
0.335 
-1.175 
-1.180 
> -1.185 

!U 

pq -1.190 
-1.195 
- 1.200 

-0.10 -0.05 0.00 0.05 0.10 

k (nm“i) 

FIG. 2. Band structure at a fixed magnetic field B = 10 Tesla 
for conduction (red) and valence band (black). The dashed 
curves show the B = 0 spectrum of Eq. (1): it is symmetric 
around the F point and consistent with Fig. 1 of Ref. 

(a). The arrows indicate possible interband transitions. 

Unlike graphene (l^ . silicene and topological 

insulators [1^, the large intrinsic band gap and the lack 
of perfect symmetry between the positive and negative 
branches of the phosphorene spectrum, shown in Fig. 2 
for B = 0 (dashed curve) and B ^ 0 (straight solid lines), 
have important implications for the peaks seen in the real 
part of cr"^(w) and the imaginary parts of tT"^(u;). In Fig. 
3 we plot the former as a function of the frequency. We 
consider a rather strong field B = 10 T so that well- 
resolved LLs are formed. The value Ep = 0.343 eV is 
between the n = 0 and n = 1 LLs, whereas the value 
Ep = 0.356 eV is between the n = 4 and n = 5 LLs. 
In the latter case the peaks for n < 3 are Pauli blocked 
and no longer possible. We notice the equally spaced 
absorption peaks. The optical selection rules allow n to 
change only by 1. In addition one needs to go from oc¬ 
cupied to unoccupied states through the absorption of a 
photon. The first peak occurring at fiuj = 1.526 eV repre¬ 
sents transitions involving the n = 0 LL. In fact, it is the 
sum of the absorption peaks of two transitions involving 


the energy differences E^ — Eq and Eq — E^ and is 
described by the last two terms in Eq. (12). This can 
also be understood from the spectrum shown in Fig. 2. 
A similar explanation holds for the peaks at fku = 1.532 
eV, fiUJ = 1.537 eV, etc. The peak spacing is propor¬ 
tional to B and experimentally one should observe such 
well-spaced peaks even for modest fields B. In contrast 
to phosphorene, in graphene and other 2D systems, the 
spectral weight of the interband peaks is redistributed 
into intraband peaks [13, [13| • This shows how the con¬ 
ductivity changes as Ep moves through the LLs. 



fioj (eV) 

FIG. 3. Real part of the longitudinal optical conductivity as 
a function of the photon energy for a field B = 10 Tesla. 



1.50 1.52 1.54 1.56 1.58 1.60 
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FIG. 4. Imaginary part of the optical Hall conductivity versus 
photon energy for a field B = 10 Tesla. 

Inspection of Eqs. (12) and (15) shows that the Hall 
and longitudinal conductivities are different due to the 
factors A®® and A®^ that reflect the difference between 
Eqs. (8) and (9). In Fig. 4 one first sees a positive 
peak and then a decrease (dip). This represents interhand 
transitions involving the n = 0 LL similar to Fig. 3 and 
the energy differences E^ — Eq and Eq — Ei[. It can 
be understood as the sum of the last two terms in Eq. 
(15). For the pure Dirac case these two peaks would 
occur at the same energy and hence would cancel out 
perfectly due to the symmetry of the spectrum. Only 
the first peak would remain in the Hall conductivity and 
all higher peaks would cancel out [ll, [ 23 | ■ 

In Fig. 5 we show the real part of the longitudinal con¬ 
ductivity (top panel) and the imaginary one of the Hall 
conductivity (lower panel) versus the field B for Huj =1.6 
eV. Due to the interband transitions from occupied to un- 
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FIG. 5. Real part of the longitudinal optical conductivity (top 
panel) and imaginary part of the optical Hall conductivity 
(lower panel) versus field B for photon energy fej = 1.6 eV. 



fico (eV) 

FIG. 6. Real part of the right-polarized optical conductiv¬ 
ity iJ+{ui) and of the left-polarized one cr_(cu) versus photon 
energy for Ef = 0.343 eV and field B — IQ Tesla. 


occupied states through the absorption of the photon, the 
oscillation patterns are similar to those shown in Figs. 3 
and 4. 

The peak structure just described above for 
and (lo) importantly affects their behaviour for 

right (-I-) and left (—) polarized light. For real experi¬ 
ments that probe the (circular) polarization of resonant 
light, as in the case of the Kerr and Faraday effects, one 
evaluates the quantity cr±(a;) given by 

u±(a;) = 5Ju-,"(a;)Tlm<(a;), (16) 

with the upp er ( lower) signs corresponding to right (left) 
polarization [I^, In Fig- 6 we show cr_(w) (solid 
black curve) and cr+{io) (dotted red curve) as functions 
of the frequency using the parameters of Figs. 3 and 4. 
As seen, there is a direct correspondence between these 
results and those of Figs. 3 and 4. The peaks in (J+{oj) 


are shifted a bit (downward) in energy relative to those 
in a-{u}). This difference also shows up in the power 
absorption spectrum given by 


P{uj) = {E/2) axx{Lo)+ayy{uj)-iayx{uj)+iaa;y{io) . (17) 


We remind that + a/// = a/// since the com¬ 
ponent = a;, y, vanishes. The component ayy{ijj) 

is given by a'/.'/.{oj) with A®’* replaced by , and 
Im tT"y(w) = —ImCTyj?(w). The spectrum P{oj) is shown 
in Fig. 7 as a function of the photon energy for two val¬ 
ues oi Ep. Given that Im(T"y(a;) is much smaller than 
5iCT"^(w), cf. Figs. 3 and 4, the peaks in it are essentially 
the same as those in the longitudinal conductivity. The 
absence of the n < 3 peaks for Ep = 0.356 is due to Pauli 
blocking and consistent with Figs. 3 and 4. 



fiaj (eV) 

FIG. 7. Power spectrum vs photon energy for an electric field 
E = 8 V/nm, for two values of i?F and field B = 10 Tesla. 


Now we consider intraband transitions between the nth 
and (n -I- l)th LLs, for Ep > 0, which involve an energy 
range much smaller than Ep. This involves large val¬ 
ues of n and is known as the semiclassical limit of the 
magneto-optical conductivity in which Ep is much larger 
than huic- Let us assume that Ep fa E/^ lies between the 
nth and (n -I- l)th LLs. The pertinent energy difference 
is Ei/i — E//_^_i = —fkjj/. For such transitions we obtain 




nd 


(w) = 7rcro^(r 


rf+ _ f+i 

1) Un+l ^(^e _ 


(18) 

The real part of cr”;)?(w) is shown in Fig. 8 by the up¬ 
per two curves. As seen, the optical spectral weight un¬ 
der these curves increases with Ep. These peaks lie in 
the range of microwave-to-THz frequencies and their am¬ 
plitude is larger than that of the interband transitions 
shown in Fig. 3. This is consistent with graphene or 
topological insulators and other symmetric 2D systems 
in which the relevant spectral weight increases with Ep, 
see, e.g., Fig. 7 of Ref. [1^, and the optical features ap¬ 
pear in the THz regime only [2^ [2^, [23 ■ The two lowest 
curves show the imaginary part of the Hall conductivity 
a'!/y{uj) for the same values of Ep. As seen, apart from 
the scale, cF//,y{u)) shows the same behaviour as cr"^(a;). 
The magnitudes of cr^y(uj) and cr^^(uj) are different due 
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ficj (meV) 


FIG. 8. Intraband limit of the real part (npper two cnrves) 
of the longitudinal optical conductivity and of the imaginary 
part (lower two curves) of the Hall conductivity versus photon 
energy for two values of Ef and B = IQ Tesla. The energy 
hw is measured from the bottom of the conduction band. 


longitudinal conductivity is the sum of the two transi¬ 
tions E~ —>• and i?~+i ^n- This corresponds 

to the first two peaks in the optical Hall conductivity 
with the first being positive and the second negative. 
The corresponding energy absorption for the first peak 
is E^ — Eq and for the second one Eq — E^. Simi¬ 
larly other peaks follow the well-defined selection rules 
between higher LLs. The spacing between the peaks also 
depends on the broadening, which we have fixed. When 
the Fermi level is in the band and not zero, the peak 
heights are suppressed and shifted downward. 


VI. SUMMARY 


to the different values of the velocity components along 
the X and y axes. The peaks of all curves occur at w = uJc 


V. OSCILLATOR STRENGTH 

The oscillator strength depends strongly on the sym¬ 
metries of the initial and final state wave functions and 
is a function of the in-plane momentum vector [s^ . The 
strength of an optical transition is typically characterized 
by the dipole matrix element between the initial and final 
states. Since a dimensionless quantity is more useful for 
making comparisons in different systems, the oscillator 
strength is often used instead of the dipole matrix ele¬ 
ment. It is defined through the A • p term in the Hamil¬ 
tonian describing the interaction between the electron 
and the electromagnetic field as 

On',n = ‘^'^{E^n'-En) |(n',s',fcy| r |n,s,fcy)|^ , (19) 

s,s' 


where m = is the mass of the particles and 

{E^, — E^)/h the frequency involved in the transitions 
from the initial to the final states. With the help of Eq. 
(3), the relevant matrix element of the position operator 
((n'l r \n) = in', s,ky\x |n, s, ky)) is 


We studied magneto-optical transport properties of 
monolayer phosphorene subject to an external perpen¬ 
dicular magnetic field. The relevant conductivities ex¬ 
hibit periodic oscillations that can be controlled by the 
magnetic field B. In each band the oscillation peaks are 
equidistant, reflecting the equally spaced LLs, and show 
a linear dependence on B. The intraband and interband 
optical transitions pertain to two completely different 
regimes: the former occur in the microwave-to-terahertz 
range and the latter ones in the visible frequency range. 
This is in contrast with a conventional 2D electron gas, 
topological insulators, and graphene in which these fea¬ 
tures appear only in the THz regime. It is also in contrast 
with phosphorene’s responses at B=0 which occur in the 
mid- to near-infrared regime 0(a), [l^, [l^. These find¬ 
ings expand the horizon of the optical properties of 2D 
phosphorene and are expected to be useful in the design 
of new optical devices. 


Appendix A 

The velocity operator, obtained from Eq. (1), reads 


dH 

9p 


a'Px + PPy 
0 


0 

->^'Px - VPy 


(A.l) 


{n,s,ky\x\n,s,ky)= xoSn'^nE 

+ {1/V2C) {V^rTTSn',n+l “b l) j ^k'y^ky- (^0) 


For s = s' the explicit evaluation of the velocity matrix 
elements gives Eqs. (8) and (9) with s = + = e and 
s = — = h. For s ^ s' we obtain explicitly 


Equations (19) and (20) clearly show the intraband 
and interband transitions shown in Figs. 2-8. The tran¬ 
sitions follow the selection rule n' = n ± 1. It is inter¬ 
esting to physically interpret the oscillator strength in 
terms of intraband and interband transitions. The intra¬ 
band transitions between the n and n -|- I states by the 
absorption of a photon shown in Fig. 8 are the same as 
those shown by Eqs. (19) and (20). For these transitions 
we have — E+ = fej®, s = s' = -b, whereas the in¬ 

terband transitions of Figs. 3-7 follow the rule n' = n± 1 
but with s ^ s'. The first absorption peak in the optical 


' ,n,n> =HVx-V^)6ky,k'y J dx4)n{u’') 

— CO 

Vn' -b I(()„'+i(m®) - Vn' (jin' -liu") . (A.2) 


Because = I.lit® we set = u® in order to have 
simple expressions. We then have 


1= Jdx(j)l^{u^)(j)niu'') PS 6n,m- (A.3) 
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This gives and „ with 


'^a:’n,n' = y/nSn,n' + l “ V^Sn,n'-l 4„.fc;,(A.4) 

'^v'^ n = ’^y} '/^Sn',n+l+VnSn',u-l 4j,.fe'-(A.5) 


To check the approximation (A.3) we evaluate 
explicitly the inte gral / for n = m {I vanishes for 
n ^ m, see Ref. [^) and n = 0,1,2,5,10 using the 
explicit expressions of the Hermite polynomials, e.g., 
Ho{x) = = 2 x,H2{x) = — 2, etc. The 

values we obtain for n = 0,1,2,5,10 are, respectively. 


0.997, 0.992, 0.983,0.924, and 0.742. This shows that the 
approximation (A.3) is a valid one at least when the 
magnetic field is strong and only a few LLs are occupied. 
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